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ABSTRACT 



In this work, the robustness of a simple, Bartlett-type processor based on matching 
broadband signal autocorrelation functions is investigated. Measures of robustness to be 
examined include the size of the localization footprint on the ambiguity surface and the 
peak-to-sidelobe levels in the presence of environmental mismatch and noise. A full-wave 
PE model is used to produce broadband replicas. Both model-generated synthetic signals, 
which provide baseline results, and measured pulses in a shallow water environment are 
analyzed. 

This work suggests that environmental mismatch has a more significant effect on 
the localization performance than noise. It also suggests that, as long as the noise level is not 
higher than the signal level, the localization performance will not be significantly affected. 
This is to be expected, since for white noise the majority of the influence on the 
autocorrelation function occurs at zero lag which has been removed in the localization 
algorithms. It is also shown that the autocorrelation matching in the time-domain is 
generally more useful for smaller bandwidths at low frequencies, which has been observed 
in previous work, whereas the autocorrelation matching in the frequency-domain is better 
suited for larger bandwidths and higher frequencies. 
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I. INTRODUCTION 



Acoustic source localization has been an intensive research issue for the past few 
decades. Prior to that, passive SONAR was used to obtain an estimate of the source 
direction and other techniques, such as Target Motion Analysis (TMA), and Contact Motion 
Analysis (CMA), were used to get an estimate of the range to the source. With the 
introduction of faster computers other possibilities arose. 

One of the techniques that was developed is known as Matched Field Processing 
(MFP). The MFP process consists of systematically placing a synthetic source at each point 
in a search grid and comparing the signal parameters estimated for all the synthetic source 
locations with the signal parameters of the true source (Tolstoy, 1993). When the synthetic 
source location and the true source location match, the correlation of the true and the 
synthetic received signal should be a maximum. Most of the work in this field has been 
done with array receivers and narrow band signals. This work deals with a single receiver 
hydrophone and broadband signals. Pioneering work in this field has mainly been done by 
Clay in the 1960’s and 70’s (Clay, 1987). 

The objective of this study is an examination of the influence of the propagation 
mismatch due to environmental uncertainty on a number of MFP algorithms for a single 
hydrophone receiver and a transient-like point source. Several aspects of environmental 
mismatch can be independently evaluated providing information of localization 
performance as a function of model degradation. Time-domain signal autocorrelation 
matching and frequency-domain autocorrelation matching are performed to produce the 
ambiguity surface for the localization. 

Measures of robustness of the algorithms to be examined include the size of the 
footprint on the ambiguity surface and the peak-to-sidelobe levels in the presence of 
environmental mismatch. A full-wave, parabolic equation model is employed to produce 
broadband replica signals. Both model -generated synthetic signals, which provide baseline 
results, and measured pulses in a shallow water environment are analyzed. A secondary but 
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no less important objective is to suggest directions for future research beyond the present 
available algorithms. 

The remainder of this thesis consists of five chapters. Chapter II describes the acoustic 
propagation model and arrival structure synthesis. Chapter El describes the matching 
algorithms and the interface of the arrival structure synthesis with these algorithms. Chapter 
IV describes the setup of the experiments and the numerical results of the matching 
algorithms. Chapter V describes the analysis of real broadband data with the matching 
algorithms. Chapter VI presents the conclusions of the study and describes issues which 
were encountered during the research and were not investigated, but which may lead to 
more robust algorithms. 
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II. THE PE-PROPAGATION MODEL AND ARRIVAL STRUCTURE SYNTHESIS 



A. THE PARABOLIC EQUATION MODEL 

In order to predict the arrival structure at a receiver location due to a point source at 
a test location, an acoustic propagation model must be employed. Throughout this work, the 
Monterey-Miami Parabolic Equation (MMPE) model (Smith, 1996) has been used to 
produce the synthetic signals. In this chapter, the general theory of the parabolic equation 
model is presented. The methods for processing the outputs of the model to produce the 
time-domain arrival structure predictions are also described. 

We begin with the inhomogeneous wave equation for the acoustic pressure p(x,t ) 
in a medium with sound speed c(3c) and density p( 3c) (Jensen et al., 1994) 



The parabolic equation model is based on an approximation to the Helmholtz wave 
equation in a cylindrical coordinate system. The majority of the ocean environment is well 
suited for a description in cylindrical coordinates. Assuming a time harmonic solution, the 
Helmholtz equation in cylindrical coordinates takes the form (see, e.g., Jensen et al., 1994) 




( 1 ) 



r dr K 



1 3 



r 



d p f (r,z,<£f \ 1 d 1 p f (r,z,<f] <? 2 p f (r,z,j) 

dr j r 2 dfi 1 dz 1 




( 2 ) 



= -4 xP 0 S(r-r s ) 



where 



p(r,z,<j>,f ,t) = p f ( r,z,<j>)e , k 0 = — is a reference wavenumber and the acoustic index 



of refraction is defined by 



n 




(3) 



3 



Note that we have assumed the density is constant in this derivation. The influence of 
density contrasts at the water-bottom interface can easily be included by defining an 
effective index of refraction which contains the appropriate additional terms (see, e.g., 
Tappert, 1977). 

As energy is primarily radiating outward from the source, p f (r, z, <t>) can be 
approximated by 

P f ( r ^,</>)= ip f (r,z,(/>)H^(k 0 r), (4) 

where vp y (r, z, 4>) is a slowly varying envelope and H^ik^ is the Hankel function of the 

first kind. In the far-field, we can take advantage of the asymptotic approximation of the 
Hankel function and re-write Eq.(4) as 

p f (r,z,( \>) = ^j=-\\i(r,z,<\>)e‘ k ° r , (5) 

normalized such that at r = R 0 ,\p f \ = P 0 . Substituting this into Eq. (2) yields 



&V , dip \d 2 ip d 2 ip 



dr 



dr r L df dz 1 



kt(n 2 -\) 



4r 



y/ = 0. 



( 6 ) 



where the source function on the right-hand side has been dropped for simplicity. The 
influence of the source term is only significant at range r = 0, and can be accounted for with 
the proper definition of the starting field. Neglecting the azimuthal coupling and the near- 
field terms, and assuming that \\i is a slowly varying function with range, we arrive at 



3\y 

dr 
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Defining the operators 



op 
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dz 2 



( 8 ) 



and 
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Eq.(7) can be written as 






( 10 ) 



With the operators defined above, this constitutes what is commonly referred to as 
the "standard" parabolic equation (SPE) (Tappert, 1977). For this work, we have employed 
the higher order "wide angle" parabolic equation (WAPE) (Thomson and Chapman, 1983) 
with operators defined by 



T = 

I WAPE 



_Lii 

ko dz 



1 + 



1 d 2 
k 2 dz\ 



Yi 



i-i 



+ 1 



(ID 



and 

U WAPE = -(n-l). (12) 

A split-step Fourier algorithm is used to numerically integrate the solution in range. This 
involves alternatively applying the U and the T op operators in the z-domain and the k z - 



domain, respectively, where each operator is simply a scalar multiplier. The algorithm for 
stepping in range from r to r + Ar can then be succinctly expressed as 

y/{r + A r,z) = e ~ ik ° ArU ^ (r ' z) ■ ^ [fFT(^*(>-,z))]* j, (13) 

where the wide angle T WAPE operator in the k 7 -space is defined as 




(14) 



Note that the FFT convention used here is consistent with the form defined by [Press, et 
al.,1988]. The calculation of the inverse FFT is accomplished by a forward FFT and a 
double conjugation. 

The output of the model is in the form of the field functions i) / (magnitude and 
phase) and has been referenced to a unity magnitude at 1 m. The field is only computed at 
discrete points and the spacing in depth and range are the primary parameters that determine 
the accuracy of the results. The wavelength of interest in this study is between 0.7 and 3 m 
and the grid spacing used is 40 cm in depth and 50 cm in range. Note that the grid size is 
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T / 

about which is necessary to obtain the highest accuracy up to ± 90 degrees 

propagation angles. This oversampling of the acoustic field was chosen to provide 
maximum information in the generation of the ambiguity surfaces. 



B. 



PE-ARRIVAL STRUCTURES 



The pressure is defined in terms of the PE field function, as stated in Eq. (6), by 

p{r,zj) = ^jrR 0 e ,k ° r \\>(r,z,f), ( 15 ) 

where P 0 is the amplitude of the source pressure measured at Rq = 1 m . The broadband 
results were obtained by running the model multiple times for all discrete frequencies in the 
bandwidth under consideration and then performing a Fourier synthesis to yield travel time 
results. In order to make this work a multiplication by some window-function is necessary. 
Assuming a normalized source amplitude, the complex arrival structure of the pressure field 
can then be written as 

p(r,z,t) = F„p{r,z,f)e-^df 

1 , ( 16 ) 

= y k ^{r,Z,f)e-^df 

i2nf— 

Since e ° r -e c ° , this phase factor can be neglected by defining the reduced time 



T = t- 



r ' 



v<V 



such that 



p(r,z,T) = ["\\i(r,z,b)e- i2,tfr df . 

Note that using a reduced time has no influence on the autocorrelation function. 



( 17 ) 
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III. LOCALIZATION ALGORITHMS 



The Transient Localization Project at the Naval Postgraduate School has in recent 
years studied different localization algorithms for the scenario of one receiver hydrophone 
and a point source. All the routines were based on fundamental concepts of generalized 
correlation functions described in Bendat and Piersol, 1971, and most of them are described 
in Miller et al., 1996, Hager, 1997, and Pierce, 1996. 

Localization algorithms may be considered generalized beamformers in which the 
plane wave replicas have been replaced by more complicated replicas of the acoustic 
propagation (e.g., modes, beams, or the vertical pressure field). The algorithms, usually 
referred to as processors, are in most cases based on a Hermitian quadratic product. The 
exact form is determined by the constraints that are put on the processor output. 

All localization algorithms are based on matching the received signal with a replica. 
Replicas are simulated received signals of synthetic sources at a large number of gridpoints 
in a search space. When the propagation model is 100% correct, and the synthetic and real 
source functions are equal, the real and synthetic received signal should exactly match 
where the real and synthetic source positions coincide. Results of the localization are 
usually presented as an ambiguity surface. All MFP algorithms are based on this principle in 
some form. 

All the localization algorithms described below are applicable to the single 
hydrophone case but can easily be extended to an array of hydrophones. The source is 
assumed to be omnidirectional. Using the reciprocity principle reduces the amount of work 
to a manageable size, and makes localization possible in a reasonable and operationally 
feasible time. The reciprocity principle states that in an environment without time variations 
(e.g., currents) the acoustic pressure at location B from an omnidirectional source at 
location A and the acoustic pressure at A from an equivalent source at B are identical if the 
densities at A and B are the same. The replicas may, therefore, be generated from 
predictions of the acoustic propagation outward from the receiver location to all possible 
source locations in the search grid. 
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A. 



TIME-DOMAIN AUTOCORRELATION MATCHING (TACM) 



Assuming a transient arrival, which can be described by the complex pressure time 
series P{t), is detected, the autocorrelation of the complex, time-domain signal is defined 
as 



T pp (T)=F a P'(t)P(t + r)dt. 



(18) 



In terms of the frequency-domain response, this can be written as (Proakis, and Manolakis, 
1996) 

P(f)=LPb)e ilnft dt, (19) 

so Eq. (18) becomes 

T„(r)^ ■ m 

Normalization of the autocorrelation produces 



r r x Vll 

pp[T) ~mf)Uf' 



( 21 ) 



Note that the autocorrelation function defined by Eq. (21) is a complex quantity. 

For a given predicted replica signal, the normalized autocorrelation function is 
similarly defined by 



Crr (*") 



MU 



( 22 ) 



Note that this function is also complex. 

The autocorrelation matching algorithm is based upon the inner product of the two 
functions defined by Eqs. (21) and (22). With proper normalization, this is 
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( 23 ) 
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Since 



r„Cj,(r)C„(r)Jr = 

= Dr 

= 1V/ 



r.df 

Hr) 

Hf) 



Hf) 



-i2nfr 



C '.df 



Hr) 



-ilnfr 



C .df 

Hf) 



Hf) 



[jr. 



" dTe- i2 * {f - f ' )r 



, ( 24 ) 



we may rewrite Eq. (23) as 



Arp — 



d f 

[r.K/)|Vr.|p(/)|V 



■ 



(25) 



Note that this quantity is real. 

In practice, the x = 0 component of the autocorrelation function is very large and 
the function defined by Eq. (25) always generates a misleading good “match”. In order to 
focus the influence on the matching of non-zero lag values, and to remove most of the 
influence of the noise (assumed uncorrelated to the signal), the x = 0 component is 
removed before forming the autocorrelation match. This is equivalent to removing the mean 
from the frequency-domain squared amplitudes, i.e.. 



Arp - 



[rjRi/)|Vr.H/)|V 



y, • 



(26) 



where the primed quantities are defined by removing the mean of the square. 



\Hff = \Hff - 



,\Hft-M- 



(27) 



and ( ) indicates an average over all frequencies. Thus, it is possible for A RP to become 

negative. Equation (26) then defines the time-domain autocorrelation matching (TACM) 
algorithm. 

Previous investigations using the TACM algorithm (de Kooter, 1997) found that the 
footprint size was very small (~ wavelength), and small, relative travel time errors (due to 
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variety of mismatch phenomena) result in large phase errors at high frequency. Thus, 
TACM seemed useful only for very low frequency signals. However, it was observed that 
while TACM might produce poor matches at correct source locations, the time-domain 

amplitude-squared signals seemed quite similar. Therefore, matching |/?(r)| 2 with (^(z)) 2 
might prove to be more robust. 



B. FREQUENCY-DOMAIN AUTOCORRELATION MATCHING (FACM) 



Matching the time-domain squared amplitude signals is equivalent to forming the 
ambiguity function 



Arp — 



rjtf(r)| 2 |p(r)| 2 <* 



r 

J-oo 



R{t) 



dtL 



Pit) 



lX 



dt 



( 28 ) 



However, based on the previous derivations, this can be defined as a matching of the 
frequency-domain autocorrelation functions. 



£ 






$ d *L 


.M* 


uf 



( 29 ) 



where 




and C SS U) 






(30) 



Note that the Fourier-Transform convention used here is chosen in order to be consistent 
with the model calculations as described in Eq.(13) (Press, et al., 1988). 

As before, we seek to reduce the overwhelming effect of the zero-lag 
(4> = 0) component by removing the mean of the squared amplitude of the time-domain 
signal. The actual form of the FACM algorithm to be implemented is then defined by 
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where 






(31) 



A 



RP 



M')f =\R(tf -(|/?(/)| 2 } ( ,|p(/)| 2 = |P'(/)| 2 -(k(/)| 2 ),, (32) 

and ( ) indicates an average over all times. Thus, it is also possible for A RP to become 
negative. 

While Eq. (31) is formally valid, there is an inherent ambiguity in the absolute 
times. In the TACM method, only relative times were used but absolute frequency was 
important. Neglecting Doppler influences, there is no ambiguity with frequency. In the 
FACM method, however, we must search over time for an optimal match. Thus, we re- 
define our FACM function as 



A rp = max, • 



f |/?'(/ + t)| 2 |p'(/)| 2 <// 

*—oo 






Vi 



(33) 



where max, represents the maximum value of the function for any value of the search 



parameter x . Therefore, by sliding the replica signal in time, we search for the best match. 
While the introduction of this new search parameter does increase the computational load of 
the algorithm somewhat, it might also have the benefit of increasing the size of the footprint 
since it may be expected that signals measured at nearby points in space decorrelate more 
rapidly in time (when no sliding is performed) than the structural character of the overall 
arrival pattern. 
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C. LOGARITHMIC AUTOCORRELATION MATCHING (LTACM & 



LFACM) 



It has recently been shown that matching logarithms of functions may be optimal for 
certain types of signals in certain environments (Makris, 1996). In addition, it was suggested 
that matching the logarithms would naturally enhance the significance of later, weaker 
arrivals. Therefore, in addition to the TACM and FACM ambiguity functions defined 
previously, we will also investigate the structure of the ambiguity surfaces defined by 



£ 


,M C 


’**w|)( io g|< 


~p,W|) 


\dx 


[f>g|^w 


OW 


1W^(t 


>i \ r 4 



(LTACM) 



(34) 



and 
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(LFACM) . 



(35) 



Note that there is now no simple relationship between these functions and their 
corresponding transform-domain functions. 
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IV. EXPERIMENTAL SETUP AND RESULTS 



To gain insight into the performance of the autocorrelation matching algorithms, several 
numerical experiments were conducted. The experiments have primarily been performed for the 
frequency band from 500 Hz to 2300 Hz. The baseline for the results was defined as the results 
of the autocorrelation matching of a PE generated source and PE generated templates. For most 
of the experiments the source range was varied between 1000 m and 10 km, the source depth 
varied between 10 m and 100 m. Received signals downrange were numerically extracted at 
depths corresponding to two short aperture arrays spanning 19.43 m to 22.36 m and 54.79 m to 
57.71 m. All ambiguity surfaces generated in this work are based on signals received at central 
elements of these arrays, either 21 m (shallow receiver) or 56.3 m (deep receiver). The generation 
of the ambiguity surfaces was performed using Matlab routines (Mathworks Inc., 1996). For the 
environments defined below, the first arrivals at this range arrive at approximately 0.8 to 7 
seconds after the transmission and the significant part of the arrival structure is about 0.4 seconds 
in length. The signals and replicas used were 2048 points in length which for a bandwidth of 
1800 Hz corresponds to a frequency spacing of 0.878906 Hz. As the project did not focus on the 
frequency characteristics of the source function, the same Gaussian spectral shape has been 
assumed for the source signal and the templates. In all presentations, the absolute scale of the 
pressure field calculation has been neglected since this drops out of the normalized ambiguity 
functions. 

Note that some plots of ambiguity surfaces displayed in this thesis are difficult to 
interpret in terms of the actual peak locations relative to the displays using a color monitor. The 
degradation is due to the limited size and print quality of the plots. 
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A. ENVIRONMENTS 



For the synthetic experiments, seven shallow water environments, including both range- 
independent and range-dependent, were defined. Typical environmental profiles for the seven 
cases are illustrated in Figs. IV. A. 1 (a-e). The synthetic signals that were used as the original 
signals for the creation of the ambiguity surfaces were generated with varying random 
uncertainties in the environmental parameters of the different scenarios. The dashed lines in Figs. 
IV. A. 1 (b-e) indicate the corresponding standard deviations. However, when generating the 
replica signals for matching, only range-independent estimates of the average environment were 
employed. 

(a) 



SSP tore*** 1 





(b) 



Figure IV.A.l: Sound speed profiles for synthetic environments (a) Case 1 and (b) Case 2 



14 



(C) 




SSP tor 4 




(d) 



SSP lor 5.«. Jnd? 




(e) 



Figure IV.A.l (continued): Sound speed profiles for synthetic environments (c) Case 3,(d) Case 
4, and (e) Cases 5-7 from SWARM experiment (Apel et al., 1997) 
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The typical parameter values used to define the bottom were a density of Q = 1.5 

b 

g/cm and an attenuation of a^= 0.10 dB/km/Hz. The water was always assumed to have a 

density of 1.0 g /cm . In addition to the third case mentioned above, an environment with the 
same parameters but with a slightly altered sound speed profile (SSP) was generated. For the 
fourth environment, replicas have been made with an average in both the SSP and the 
bathymetry, with an average in the SSP only, with an average in the bathymetry only, and 
with all parameters known exactly. Replicas have also been made for the first environment 
with a different bottom bathymetry assuming c =1575m/s with a gradient of 1.5 (m/s)/m in 
contrast to the bottom parameters mentioned above. 

Case 1 is the only scenario where the environment was both range-independent and 
exactly known. In all other cases, environmental uncertainties were introduced. Cases 5 
through 7 have a similar environment defined by a field of nonlinear internal solitons of 
varying extent within the water column. The differences between these cases were due to the 
placement of the source and receiver within the soliton field. The variability of sound speed 
due to the presence of these solitons is indicated in Fig. IV.A.l.(e). The respective 
uncertainties introduced for the different cases will be stated in detail in the next section 
together with the discussion of the results for the different matching algorithms. 

For each of the seven scenarios above, four different source locations and two 
different receiver locations were chosen. The bottom bathymetry was treated as range- 
independent for the generation of all replica templates although some minor bathymetry 
fluctuations were used in the generation of the original synthetic signal for most cases. 



B.l. Case 1: Isospeed sound speed profile, 100 m deep 

In this case the environment was exactly the same for the generation of the replica 
signal and the synthetically generated measured signal. Synthetic data for four different source 



16 



locations between 2.6 km and 9.9 km in range and 18 m to 83 m in depth were investigated 
at receiver depths of 21 m and 56.3 m. An analysis of the correlation between propagation 
loss and pulse spread of the received signal with respect to the "goodness" of a localization 
did not provide any meaningful results. The maximum spatial cross-correlation of the 
frequency-domain autocorrelation functions of the signals across the shallow array as seen in 
Fig. IV.B.l(a), however, shows that for this scenario there are distinct differences between 
these functions at different depths even though the corresponding receivers were only 
seperated by about 20 cm. This high spatial decorrelation may infer a higher probability of a 
good, less ambiguous localization since each array element receives a unique arrival structure. 

For all four source locations and both receiver depths, distinct localizations could be 
made at the correct source positions. Some results are shown in Figs. IV.B.l (b-f). The 
TACM method provided a larger footprint size when only the low frequencies, i.e. the first 
quarter (500-950 Hz), of the total bandwidth (500-2300 Hz) were considered as displayed in 
Figs. IV.B.l.(b) and (c). The footprint size is measured as the area adjacent to the location of 
the peak value in which the corresponding values on the ambiguity surface equal at least half 
the peak value. When all frequencies were used, the footprint was roughly 3.5 m x 25 m as 
compared to a footprint of nearly 4.5 m x 55 m when using only the first quarter of the 
bandwidth. This is consistent with previous work (de Kooter, 1997) which showed that 
TACM is highly sensitive to phase mismatch at high frequencies. 

For the FACM method, the use of only the lower or higher frequencies of the total 
bandwidth affected the results negatively providing a smaller footprint size and higher 
sidelobe levels. The FACM method produced basically identical results when only every 
fourth data point was sampled of the full bandwidth. This is to be expected, however, since 
the frequency sampling was roughly four times smaller than needed to avoid aliasing or wrap- 
around in the time domain. The FACM generally provided a larger footprint size than the 
TACM, roughly 5 m x 80 m. This was especially obvious when the search grid was 
drastically undersampled spatially. The best spatial sampling would be provided by using a 
grid spacing of about one half of the minimum wavelength present. For a sound speed of 
1500 m/s and a maximum frequency of 2300 Hz, we find the best spatial sampling at 0.3261 
m. Thus all scenarios that were analyzed were spatially undersampled, some extremely but 
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others only slightly. 

The FACM always provided the larger peak values as compared to the TACM. The 
peak-to-sidelobe levels, however, were high (> 4 dB) for all localizations for both algorithms 
even when the environment was spatially undersampled. The location of the receiver did not 
have any significant effect on the results. They were almost identical for both shallow and 
deep receiver locations. 

Changing the sound speed in the bottom to a higher value and using a different bottom 
sound speed gradient for the generation of the measured signals did almost not affect the 
results at all. For all four source locations and both receiver depths, distinct localizations 
could be made at correct source positions using both algorithms with the same footprint-size 
and peak-to-sidelobe levels as before. 



(a) 




Figure IV. B. 1: (a) Maximum spatial correlation of the frequency-domain autocorrelation 
functions across the shallow array at 5.3 km. 
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(b) 
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Figure IV.B.l (continued): TACM for Case 1. (b) Ambiguity surface for a shallow receiver 
localizing to a source location at 5.3 km, 33 m deep. Spatial sampling 128 points in depth and 150 
points in range. All 2048 data points are used.(c) Ambiguity surface for shallow receiver 
localizing to a source location at 5.3 km, 33 m deep. Spatial sampling 128X150. Only first 512 
data points are used (bandwidth 500-950Hz). 
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(d) 





Figure IV.B. 1 (continued): TACM for Case 1. (d) Ambiguity surface for deep receiver localizing 
to a source location at 7.3 km, 61 m deep. Spatial sampling 32X300. All 2048 data points are 
used. 




Figure IV.B. 1 (continued): FACM for Case l-(e) Ambiguity surface for shallow receiver 
localizing to a source location at 5.3 km, 33 m deep. Spatial sampling 128X150. 512 data points 
are used (full bandwidth). 
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Figure IV. B. 1 (continued): FACM for Case 1. (f) Ambiguity surface for shallow receiver 
localizing to a source location at 5.3km, 33 m deep with mismatch in bottom parameters. Spatial 
sampling 32X300. 512 data points are used (full bandwidth). 



B.2. Case 2: Isospeed sound speed profile, 200 m deep 

As in Case 1 four different source locations with both receiver depths were analyzed. In 
this case however, a slight environmental mismatch was introduced. Instead of a strictly isospeed 
SSP, an uncertainty of (+/- 3 m/s) in the SSP was applied to a few depth points in the SSP for the 
generation of the synthetic measured signals. The replicas were then generated using only an 
average or approximate SSP. Correct localizations could be made in all cases using both 
algorithms and both receiver depths even when the environment was spatially undersampled. The 
results are comparable to those found for Case 1. However, a slight decrease in the peak-to- 
sidelobe levels was observed compared to Case l, more so for the TACM (3 dB) than for the 
FACM (4 dB). The maximum values also decreased as compared to Case 1 which was expected 
because of the environmental uncertainty. The footprint size for the FACM remained almost 
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unchanged compared to Case 1 whereas it decreased for the TACM to approximately 3.5 m x 45 
m. Again, the FACM proved to be more robust than the TACM to spatial undersampling and 
environmental mismatch. As in Case 1 the location of the receiver did not significantly affect the 
results. The spatial correlation of the frequency domain autocorrelation functions was 
comparable to Case 1 . 



B.3. Case 3: Positive sound speed gradient 

Again four different source locations with both receiver depths were analyzed. In this 
case, an uncertainty in the bottom bathymetry was introduced. Specifically, the depth was varied 
by (+/- 3 m) with range from the original 100 m for the generation of the measured signal. A 
further decrease in the peak-to-sidelobe levels was observed for both the FACM and TACM as 
compared to Case 2. The footprint size for the FACM did not change much compared to Case 2, 
as seen in Fig. IV.B.3(d), while it decreased further for the TACM to 3 m x 35 m, as displayed in 
Fig. IV.B.3(b). Correct localizations, however, could still be made in all cases with both 
algorithms despite undersampling spatially. 

When a mismatch in the sound speed profile was introduced instead of the uncertainty in 
the bathymetry, the results did not change significantly with respect to the shallow receiver. 
However, the localization with the deep receiver became much poorer for both algorithms. The 
peak-to-sidelobe levels decreased drastically to approximately 1 dB for the TACM and 2 dB for 
the FACM while the footprint size decreased to 3 m x 30 m for the TACM and 6 m x 64 m for 
the FACM as seen in Figs. IV.B.3(c) and (e). The worst results occured when both source and 
receiver were at a deep location. This can be explained by the fact that in this environment the 
acoustic rays are bent upwards so that at a distant deep receiver the pressure field is more 
ambiguous due to fewer arrivals with respect to the possible source location. Localizations, 
however, were still made even when the search grid was spatially undersampled. 

For both algorithms, the mismatch in the SSP had a stronger effect on the results than the 
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mismatch introduced in the bathymetry. Overall the FACM proved to be more consistent for the 
same source locations than the TACM. It can be stated that an environmental mismatch in more 
complex SSP scenarios (as in Case 3) affects the results more than the same uncertainty in a less 
complex scenario (as in Case 2). 

The maximum spatial cross-correlation of the frequency-domain autocorrelation 
functions across the shallow array, as seen in Fig. IV. B. 3(a), illustrates a slight increase in the 
correlation of the signals as compared to Case 1 and Case 2. This is consistent with the observed 
increasing ambiguity (i.e., lower peak-to-sidelobe levels) for the localization process. Thus, the 
higher the spatial correlation, the more difficult it becomes to determine the true source 
localization. 




Figure IV. B. 3: (a) Spatial correlation of the frequency-domain autocorrelation functions across 
the shallow array at 7.2 km. 
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Figure IV.B.3 (continued): TACM with mismatch in only the bathymetry, (b) Ambiguity surface 
for a deep receiver localizing to a source location at 9.3 km, 59 m deep. Spatial sampling 32X300. 
All 2048 data points are used. 



(c) 





Figure IV.B.3 (continued): TACM for Case 3 with mismatch in sound speed profile only, (c) 
Ambiguity surface for a deep receiver localizing to a source location at 7.2 km, 53 m deep. Spatial 
sampling 32X300. AH 2048 data points are used. 
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Figure IV.B.3 (continued): FACM for Case 3 with mismatch in bathymetry only.(d) Ambiguity 
surface for deep receiver localizing to a source location at 9.3 km, 59 m deep. Spatial sampling 
32X300. Only 512 data points are used (full bandwidth). 



(e) 




Figure IV.B.3 (continued): FACM for Case 3 with mismatch in SSP only, (e) Ambiguity surface 
for deep receiver localizing to a source location at 7.2 km, 53 m deep. Spatial sampling 32X300. 
512 data points are used (full bandwidth). 
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B.4. Case 4: Strong negative gradient in SSP 



Again four different source locations were analyzed for both receiver depths and both 
algorithms, FACM and TACM. In additon to a fairly high uncertainty in the sound speed profile, 
an uncertainty in the bathymetry was introduced. For the TACM, correct localizations could only 
be made when the search grid was spatially well sampled. Even then the peak had to be enhanced 
by raising all values for the ambiguity surface to a higher power to make it distinctively visible in 
the plot (equivalent to re-defining the color palette to enhance lower values) as seen in Fig. 
IV.B.4(b). There was no significant improvement observed in the results for the TACM when 
considering only the low frequency band or the full bandwidth due to the very strong 
environmental mismatch. 

For the FACM, localizations could be made in most cases despite a rather poor footprint 
size of 3 m x 32 m and low peak-to-sidelobe levels around 1 dB as seen in Fig. IV.B.4(d). This 
was true also for the strongly spatially undersampled search grid. Sampling only every fourth 
point of the given bandwidth produced the same results as compared to sampling all points, as 
before. 

A technique that seemed to provide improved results for the FACM was zero-padding the 
signal in the frequency-domain to twice its length. By doing this, the time-domain spacing is 
reduced to half and therefore additional data points, i.e. a better time resolution, is provided 
(Ziomek, 1995). Thus, the probability of a correct localization by sliding the time window over 
twice as many data points increases. Of course, this provides no true additional information. But 
since increments in time are closely related to increments in range for such propagation 
problems, the influence of this zero-padding technique is similar to an increased spatial 
resolution which the FACM algorithm takes advantage of by searching for the maximum match 
when sliding the time-domain signal. Additional zero-padding would not be expected to 
continually improve the results, and so was not considered further. Also, the computational effort 
increases when zero-padding a signal before matching it with the replicas. 

Isolating the uncertainties, i.e. considering a mismatch in the sound speed profile only 
with no mismatch in bathymetry and then vice versa, showed for both algorithms that a mismatch 
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in the sound speed profile is much more devastating to the localization results than a mismatch in 
the bathymetry. The worst results occured when both types of mismatch were introduced. In most 
cases, the FACM proved again to be more robust and capable of providing better results than the 
TACM as displayed in Figs. IV.B.4(b) and (d). However, for many cases with a mismatch in only 
the SSP, no localization could be made using either algorithm. The ambiguity surface showed too 
many false matches with values close to one, so that a distinction between peak and sidelobes 
could no longer successfully or correctly be made. For both algorithms, the results improved 
considerably with improved spatial sampling as seen in Figs. IV.B.4(c),(e). It is obvious that 
when the environmental mismatch increases the spatial sampling becomes more important and 
needs to be increased in order to provide good results. 

Generally, the very strong negative gradient in the soundspeed profile and the resulting 
strong downward bending of the acoustic energy makes a localization extremely difficult. 

The maximum spatial cross-correlation of the frequency-domain autocorrelation functions across 
the deep array, as seen in Fig. IV.B.4(a), shows a rather poor decorrelation (relatively high 
correlation) between these quantities at the corresponding receiver locations, which infers an 
even more difficult, i.e. more ambiguous, localization process as compared to Case 3. 



(a) 




Figure IV.B.4 (a): Spatial cross-correlation of the frequency-domain autocorrelation functions 
across the deep array at 7.2 km. 
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Figure IV. B. 4 (continued): TACM for Case 4. (b) Ambiguity surface for shallow receiver 
localizing to a source location at 7.2 km, 18 m deep. Spatial sampling 32X300. All 2048 data 
points are used. No match. 




Figure IV.B.4 (continued): TACM for Case 4 with mismatch in sound speed profile and 
bathymetry, (c) Ambiguity surface for shallow receiver localizing to a source location at 5.7 km, 
31 m deep. Spatial sampling 128X150. All 2048 data points are used. 
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(d) 




Figure IV.B.4 (continued): FACM for Case 4 with mismatch in SSP and bathymetry, (d) 
Ambiguity surface for shallow receiver localizing to a source location at 7.2 km, 18 m deep. 
Spatial sampling 32X300. 512 data-points used. Values raised to the fifth power. 




(e) 



Figure IV.B.4 (continued): FACM or Case 4. (e) Ambiguity surface for shallow receiver 
localizing to a source location at 5.7 km, 31 m deep. Spatial sampling 128X150. Every fourth 
data point is used. Mismatch in SSP and bathymetry. Values raised to fifth power. 
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B.5. Case 5: Realistic, slightly range-dependent environment 



For the generation of the synthetic measured signals in this case, a moderately range- 
dependent bathymetry similar to that of Case 3 was used. In addition, the sound speed structure in 
the water column was extracted from a real experimental database provided by Chiu (Chiu et al., 
1997). This data was taken from a shallow water experiment and includes the influence of a 
nonlinear internal soliton wave packet which can significantly perturb the sound speed structure, 
as seen in Fig. IV. B. 5(a). The sound speed information (in m/s) is given by the color coding. For 
this case, the sound speed profiles were extracted between 25 km and 35 km of this database. The 
result is a slightly range-dependent SSP with no major perturbations due to the soliton packets. 
For the generation of the synthetic data sets, the position in range of the receiving arrays was 
assumed fixed while the various source locations were changed in range and depth. For the 
generation of the replicas, only an average bottom depth and average SSP were used to produce a 
range-independent environment. Thus, mismatch existed in bathymetry and sound speed, and the 
"real" environment was range-dependent in both while the replicas assume an approximate range- 
independent scenario. 

Despite these uncertainties and a rather poor spatial decorrelation of the frequency-domain 
autocorrelation functions at the receivers as displayed in Fig. IV.B.5(b), the results for both 
algorithms, TACM and FACM, were similar as seen in Figs. IV.B.5(c)-(f). However, the relative 
position of source and receiver proved to be important for the localization. Specifically, due to the 
presence of a sound channel, the similar location of either both source and deep receiver in the 
channel or both source and shallow receiver above the channel provided good localizations. On 
the other hand, in cases where only either source or receiver were put into the SOFAR channel 
the results were poor. No significant differences between the localizations made using FACM and 
TACM were observed apart from a better footprint size for the FACM (6 m x 64 m) as compared 
to the TACM (3 m x 32 m). The peak-to-sidelobe levels were approximately 2 dB for both 
algorithms. 
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Figure IV. B. 5; (a) Soliton field between approximately 5 and 30 m depth and 5 and 20 km range 



(b) 




Figure IV. B 5 (continued): (b) Spatial cross-correlation of the frequency-domain autocorrelation 
functions across the deep array at 4.4 km. 
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(d) 



Figure IV.B.5 (continued): TACM for Case 5 with mismatch in sound speed profile (range 
dependence), (c) Ambiguity surface for deep receiver localizing to a source location at 5.9 km, 5 1 
m deep. Spatial sampling 171X500. All 2048 data points are used, (d) Ambiguity surface for a 
deep receiver localizing to a source location at 7.2km, 32m deep (just within sound channel). 
Spatial sampling 171X500. All 2048 data points are used. 
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Figure IV.B.5 (continued): FACM for Case 5. (e) Ambiguity surface for deep receiver localizing 
to a source location at 5.9 km, 51 m deep. Spatial sampling 171X500. 512 data points are used, 
(f) Ambiguity surface for deep receiver localizing to a for source location at 7.2 km, 32 m deep. 
Spatial sampling 171X500. Every fourth data point is used (full bandwidth). 
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B.6. Case 6: Range-dependent SSP pertubed by internal wave soliton 
field from 7.5 km outward 



For the generation of the synthetic signals in this case, the SSP was made much more 
range-dependent by extracting from the experimental data set between ranges 0 km and 10 
km (refer to Fig. IV.B.5(a)). This produced an environment perturbed by the soliton packet 
from about 7.5 km outward from the receiver array. As before, the environment was assumed 
fixed relative to the receiver position while the various source locations were moved about 
such that the source may lie within the range of soliton influence. Again, a single, average 
SSP and bottom depth were used in the generation of the replicas. 

In this case, the peak-to-sidelobe levels were a little lower for both algorithms (1 dB) 
as compared to Case 5 for all source locations including those closer than the soliton field. 

For those source locations closer than the soliton field, this is presumably due to the influence 
of the solitons on the average, range-independent SSP. For the source locations within the 
range of the soliton field the localization was excellent for both source and receiver in the 
channel using both FACM and TACM. Again, the footprint-size for the FACM was larger 
than for the TACM and approximately the same as in Case 5 as displayed in Figs. 
IV.B.6(a),(c),(d),(f). For the shallow receiver, however, no localization was made for the 
TACM whereas the FACM provided a rather poor but correct localization as seen in Figs. 
IV.B.6(b) and (e). Therefore, the FACM continues to be more robust against environmental 
mismatches than the TACM, being able to make correct localizations when the TACM fails. 

The presence of the SOFAR channel represents a rather favorable SSP for good 
localizations, and thus the environmental mismatches introduced can be larger than those in 
Case 4 and still provide correct localizations. The cross-correlation between the frequency- 
domain autocorrelation functions of the receivers of the deep array are effectively the same as 
for Case 5, suggesting that little or no energy from the synthetic source locations interacted 
with the soliton field. The same comments then apply. 
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Figure IV. B, 6: TACM for Case 6 with range dependent SSP. (a) Ambiguity surface for deep 
receiver localizing to a source location at 8.8 km, 75 m deep. Spatial sampling 171X150. All 2048 
data are points used.(b) Ambiguity surface for shallow receiver localizing to a sourcelocation at 
8.8 km, 75 m deep. Spatial sampling 171X150. All 2048 data points are used. No match. 
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Figure IV. B. 6 (continued): TACM for Case 6 with range dependent SSP. (c) Ambiguity surface 
for deep receiver localizing to a source location at 2.4 km, 47 m deep. Spatial sampling 171X150. 
All 2048 data points are used. 




Figure IV. B. 6 (continued): FACM for Case 6. (d) Ambiguity surface for deep receiver localizing 
to a source location at 8.8 km, 75 m deep. Spatial sampling 171X150. All data points are used. 
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(C) 






Figure IV.B.6 (continued): FACM for Case 6. (e) Same as (d) with shallow receiver. Values 
raised to the fifth power, (f) Ambiguity surface for deep receiver localizing to a source location at 
2.4 km, 47 m deep. Spatial sampling 171X500. Every fourth data point is used. Values raised to 
the fifth power. 
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B.7. Case 7: Range-dependent SSP perturbed by internal wave soliton field 
from 2.5 to 8 km 

The sound speed structure used to produce the synthetic measured signals for this case 
was extracted from the experimental data set between ranges 5 km and 15 km (refer to Fig. 
IV.B.5(a)). This produced an environment perturbed by the soliton packet from about 2.5 km 
to 8 km from the receiver arrays. Again, the environment was fixed relative to the receiver 
range, and only an average, range- independent SSP and bottom depth were used in the 
generation of the replica fields. 

For source locations closer than the soliton field and above the SOFAR channel, the 
localizations were rather poor for both TACM and FACM as seen in Figs. IV.B.7(a) and (d). 
However, the footprint-size and peak-to-sidelobe levels were once again better for the FACM 
(4 m x 64 m, 1 dB) than for the TACM (3 m x 32 m, 0.4 dB). For the TACM method, the 
values on the ambiguity surface had to be raised to the fifth power in order to make the peak 
distinctively visible. This underscores the superiority of the FACM over the TACM in the 
presence of mismatch. 

For the source locations beyond the range of the soliton field and within the SOFAR 
channel, a good localization could be made with both FACM and TACM as displayed in 
Figs. IV.B.7(c) and (f). In this unusual case, the TACM even provided slightly better results. 

For the source locations within the soliton field, the closer source could be located, as 
seen in Figs. IV. B. 7(b) and (e), whereas the source further away could not be localized by 
either algorithm. As in Case 5 and Case 6 the relative source and receiver locations were 
important for the goodness of the match. 

Generally, it can be stated that the presence of such strongly range-dependent 
perturbations does affect the results significantly, and that the results get worse the larger the 
perturbations are. 
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Figure IV.B.7: TACM for Case 7 with mismatch in sound speed profile (range dependence), (a) 
Ambiguity surface for shallow receiver localizing to a source location at 1.8 km, 20 m deep 
(before soliton field). Spatial sampling 171X500. All 2048 data points are used. Values raised to 
fifth power, (b) Ambiguity surface for deep receiver localizing to a source location at 3.2 km, 58 
m deep (within soliton field). Spatial sampling 171X500. All 2048 data points used. 



39 






Figure IV. B. 7 (continued): TACM for Case 7 with mismatch in sound speed profile (range 
dependence), (c) Ambiguity surface for deep receiver localizing to a source location at 8. 1 km, 42 
mdeep (just behind soliton field). Spatial sampling 171X150. All 2048 data points used. 




Figure IV. B. 7 (continued): FACM for Case 7. (d) Ambiguity surface for shallow receiver 
localizing to a source location at 1.8 km, 20 m deep. Spatial sampling 171X150. All data points 
are used. Values raised to the fifth power. 
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Figure IV.B.7 (continued): FACM for Case 7. (e) Ambiguity surface for deep receiver localizing 
to a source location at 3.2 km, 58 m deep. Spatial sampling 171X500. Every fourth data point is 
used. 
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Figure IV.B.7 (continued): FACM for Case 7. (f) Ambiguity surface for deep receiver localizing 
to a source location at 8. 1 km, 42 m deep. Spatial sampling 171X500. Every fourth data point is 
used. Values raised to the fifth power. 
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C. LOGARITHMIC FREQUENCY-DOMAIN AND TIME-DOMAIN 
AUTOCORRELATION MATCHING 



In an attempt to enhance the significance of later, weaker arrivals, logarithmic versions of 
TACM and FACM were implemented (see Chapter III). However, it was found that such an 
approach degraded previous results where successful localizations had been made. Figure C. 1 
shows the results obtained for the easiest case (Case 1). It is obvious that for this case, a perfectly 
known and spatially well sampled environment, the footprint size and the peak-to-sidelobe levels 
are much worse using the logarithmic algorithms than the standard TACM and FACM For this 
case, the peaks can be enhanced so that a distinct localization can be made. But in all other cases 
no match could be found with either logarithmic algorithm. Therefore, no further results from 
these logarithmic algorithms will be discussed. 




Figure IV.C. 1 : Logarithmic Autocorrelation Matching, (a) Frequency-Domain Autocorrelation 
Matching. Ambiguity surface for shallow receiver localizing to a source location at 5.3 km, 33 m 
deep. Spatial sampling 128X150. All data points are used. 
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Figure IV.C. 1 (continued): Logarithmic Autocorrelation Matching, (b) Same as (a) for Time- 
Domain Autocorrelation Matching. 



D. Influence of noise on the performance of algorithms 



So far, all of the synthetic data discussed in this chapter have been noise-free. This section 
investigates the effect of adding white noise to the synthetic data with regard to the consistency of 
the algorithms. In other words, how much noise can be added to the signal so that a localization is 
still possible in those cases where successful localizations occurred in the absence of noise? For 
this analysis, three different cases where localizations could originally be made for both TACM 
and FACM were investigated. The level of added white noise was varied such that the signal-to- 
noise ratio (SNR) was reduced from 20 dB down to - 20 dB. The SNR is defined as the average 
signal power divided by the average noise power (Oppenheim, 1983). Figure IV.D. 1 shows the 
time-domain and frequency-domain signals with the added noise at the receiver locations for Case 
1 . 
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Figure IV.D.l: (a) Time-signal 0 to 1.38 sec. at receiver for source at 5.3 km 33 m deep for Case 
1 (isospeed SSP) for different signal-to-noise ratios, (b) Frequency spectrum (500-2300 Hz) of 
time-signal. 



For Case 1 where the environment is perfectly well known and the search grid spatially 
well sampled, the results show that for the FACM algorithm a good localization can still be made 
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for an SNR down to -3 dB. At an SNR of -10 dB, no localization can be made. For the TACM 
algorithm, good localizations can be made down to an SNR of 0 dB. From -3 dB down, the 
localization becomes very poor to impossible. This confirms once again that the FACM is more 
robust against uncertainties than the TACM as seen in Fig. IV.D.2. 





Figure IV.D.2: (a) TACM. Ambiguity surfaces for source at 5.3 km, 33 m deep for Case 1 
(isospeed SSP) for different levels of SNR. (b) FACM. Ambiguity surfaces for source at 5.3 km, 
33 m deep for Case 1 (isospeed SSP) for different levels of SNR. 
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For Case 2 with a nearly isospeed SSP and an uncertainty in the SSP, the results are 
basically the same as for Case 1. Again, the localization works well down to -3 dB SNR for the 
FACM, which performs better than the TACM, as displayed in Fig. IV.D.3. It can also be seen 
that the results in Case 2 are worse as compared to Case l for the same SNR's which is due to the 
fact that the search grid is spatially undersampled and a slight mismatch exists in the SSP. 





Figure IV.D.3: (a) TACM. Ambiguity surfaces for source at 3.1 km, 22 m deep for Case 2 
(isospeed SSP, uncertainty in SSP) for different levels of SNR. (b) FACM. Ambiguity surfaces for 
source at 3. 1 km, 22 m deep for Case 2 (isospeed SSP, uncertainty in SSP) for different levels of 
SNR. 
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The last case to which noise was added was one data set from Case 7 which contains 
a mismatch in bathymetry as well as a range-dependent mismatch in SSP due to the presence 
of a packet of solitons. The specific data set examined was one for which only a poor 
localization with low peak-to-sidelobe levels was possible for both TACM and FACM 
without noise. For the FACM a reasonable localization could still be made for -3dB SNR 
despite the significant environmental mismatches already present, as seen in Fig. IV. D. 4. All 
of these results suggest that the influence of environmental mismatch has a more significant 
effect on the localization results than the presence of noise. It also suggests that, as long as 
the noise level is not higher than the signal level, the results will not be significantly affected 
by the presence of noise. This result is not entirely unexpected since the majority of the noise 
influence in the autocorrelation function occurs at zero lag which has been removed in both 
TACM and FACM algorithms. In the real world, a detected signal in noise will usually not 
have negative SNR. Otherwise, a detection would be rather improbable in the first place. 

For the TACM, the results were similar to the FACM although continuing to show 
less robustness than the FACM, as displayed in Fig. IV.D.4. The TACM fails to make 
localizations at an SNR of -3 dB whereas the FACM still worked well. At an SNR of 0 dB, 
the TACM localization is very poor. 
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Figure IV.D.4: (a) FACM. Ambiguity surfaces for source at 1.8 km, 20 m deep for Case 7 with 
mismatch in SSP and bathymetry for different levels of SNR. 
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OdB SNR 




Figure IV.D.4 (continued): (b) TACM. Ambiguity surfaces for source at 1.8 km, 20 m deep for 
Case 7 (Mismatch in SSP and bathymetry) for different levels of SNR. 
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E. Effect of bandwidth in the presence of noise 



In this section, the effect of a reduction in the bandwidth on the results as obtained in 
Section IV. D is analyzed. The scenario of Case 2 was chosen as a representative environment. For 
both algorithms, TACM and FACM, a good localization was still possible for the spatially 
undersampled case with an SNR of -3 dB (see previous section). Thus the noise-level used in this 
experiment was fixed to -3 dB SNR, and the bandwidth was varied for both algorithms to find the 
limits of a possible localization. 

Figures IV.E. l(a-h) show the results for the TACM. The TACM worked best, as shown 
previously, for lower frequency bands. A localization was still possible for a bandwidth of only 
225 Hz, from 500 Hz to 725 Hz, which is equal to just one eighth of the full bandwidth with 256 
data points as seen in Figs. IV.E. 1(a) and (b). Considering only the lowest sixteenth of the 
bandwidth, 500 to 612.5 Hz, led to no localization. This may be due to fewer arrivals being 
adequately resolved at smaller bandwidths which causes the autocorrelation and its match to drop 
off rapidly away from the source location. 

When sampling the full bandwidth at only every fourth and every eighth point, 
localizations could still be made despite being rather poor. For sampling only every sixteenth data 
point of the full bandwidth, there was still a relative peak at the correct location but the absolute 
maximum was found elsewhere as seen in Fig. IV.E. 1(c). The undersampling introduces aliasing 
effects and severe wrap-around in the time-domain of the autocorrelation functions so that a 
localization with these degraded functions becomes impossible. Looking at only the lower half 
and the upper half of the bandwidth and sampling only every fourth data point still provided a 
correct localization, although footprints and peak-to-sidelobe levels were very poor, as seen in 
Figs. IV.E.l(d-h). 

Summarizing the results for the TACM, it can be stated that it works best in the low 
frequency region, and that a bandwidth of only a few hundred Hertz (in this case one eighth of the 
total bandwidth) is sufficient to provide reasonable localization matches even for a SNR of -3 dB. 
It should not be forgotten, however, that the environmental mismatch in this case was only of 
little significance. 
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Figure IV.E. 1: TACM. Ambiguity surfaces for source at 3. 1 km, 22m deep for Case 2 with SNR 
of - 3dB. (a) first quarter of full bandwidth (500-950Hz), (b) first eighth of full bandwidth (500- 
725Hz), (c) first sixteenth of full bandwidth (500-612.5 Hz), no match. 
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Figure IV.E. 1 (continued): TACM. Ambiguity surfaces for source at 3. 1 km, 22 m deep for Case 
2 with SNR of - 3dB. (d) lower half of full bandwidth, every fourth data point is used, (e) upper 
half of full bandwidth, every fourth data point is used. 
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Figure IV.E. 1 (continued): TACM. Ambiguity surfaces for source at 3. 1 km, 22 m deep for Case 
2 with SNR of - 3dB. (f) full bandwidth, every fourth data point is used, (g) full bandwidth, every 
eighth data point is used, (h) full bandwidth, every sixteenth data point is used. 

For the FACM, the results show that the best localizations can be made when the full 
bandwidth of the signal is exploited, as suggested previously. Again, the Case 2 scenario was 
chosen as a representative environment. When sampling down to every eighth point, good 
matches are still obtained as seen in Figs. IV.E.2(a-c). 

To examine the performance of the FACM algorithm for smaller bandwidths, calculations 
were made for half of the original bandwidth including the lower half (500 Hz- 1400 Hz), the 
middle half (950 Hz- 1850 Hz), and the upper half (1400 Hz-2300 Hz). The results, as displayed in 
Figs. E2 (d)-(f), illustrate that the FACM algorithm works much better for the higher frequencies 
than for the lower frequencies. This may be due to the increased complexity at higher frequencies 
leading to more unique structures in the autocorrelation function. The best results, however, are 
reached when sampling the middle of the bandwidth. Finally, when only one quarter of the 
bandwidth is used, the results become much worse as seen in Figs. IV.E.2(g-i). For smaller 
bandwidths, no localizations could be made. 
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Figure IV.E.2: FACM. Ambiguity surfaces for source at 3. 1 km, 22 m deep for Case 2 with SNR 
of - 3dB. (a) full bandwidth, every fourth data point is used, (b) full bandwidth, every eighth data 
point is used, (c) full bandwidth, every sixteenth data point is used. 
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Figure IV.E.2 (continued): FACM. Ambiguity surfaces for source at 3. 1 km, 22 m deep for Case 
2 for SNR of - 3dB. (d) lower half bandwidth, every second data point is used, 

(e) middle half bandwidth, every second data point is used, (f) upper half bandwidth, every second 
data point is used. 
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(h) 
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Figure IV. E. 2 (continued): FACM. Ambiguity surfaces for source at 3. 1 km, 22 m deep for Case 
2 with SNR of - 3dB. (g) upper quarter bandwidth, every second data point is used, no match, (h) 
same as (g) with all data points used, (i) middle quarter bandwidth, all data points are used. 
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F. Summary of synthetic data results 

The FACM algorithm provided a larger footprint size and higher peak-to-sidelobe levels in 
nearly all synthetic cases as compared to the TACM. The TACM is generally more useful for 
smaller bandwidths at low frequencies, whereas the FACM is better suited for larger bandwidths 
and higher frequencies. The FACM proved to be more robust against environmental mismatches, 
added white noise, and spatial undersampling than the TACM. The influence of environmental 
mismatch has a much more significant effect on the localization results than the presence of noise 
for both algorithms. This is to be expected, since for white noise the majority of the influence on 
the autocorrelation function occurs at zero lag which has been removed in both TACM and 
FACM algorithms. Good localizations even in the presence of severe environmental mismatch 
could still be made with both algorithms for SNR's down to at least -3 dB. 

For good localization results the knowledge of the SSP is more crucial than the estimate 
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of any other environmental parameter. Furthermore, the best localization results are obtained in 
isospeed environments followed by environments with a SOFAR channel. In this case the 
relative location of source and receiver are critical for good results. Good localizations could be 
made with either both source and receiver in the channel or both source and receiver outside the 
channel. Environments characterized by positive, upward refracting sound speed gradients are 
generally more favorable than environments with negative, downward refracting sound speed 
gradients which represent the most difficult localization environment. 

The FACM provided correct localizations for almost all cases and for many scenarios in 
which the TACM was unsuccessful. Improved spatial sampling or techniques such as zero- 
padding the signal to twice its length which, of course, increase the necessary computational 
time, improved the results for the FACM in the more difficult scenarios. 
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V. REAL DATA RESULTS 



In this section, the results of the analysis of some real broadband explosive SUS data 
taken from the Mid-Atlantic Bight Field Study (Pickart, et ah, 1996) are presented. The Mid- 
Atlantic Bight Field Study, also known as Primer experiment, was a collaborative shallow 
water acoustics experiment between the Woods Hole Oceanographic Institute, the University 
of Rhode Island, and the Naval Postgraduate School. The purpose was to study the influence 
of a highly variable, strongly 3-D, shallow water environment on acoustic propagation. 

Among the assets used were 61 P-3 SUS deployments, and two 16-element vertical receiver 
arrays. For this work, recorded data of a SUS detonation about 6.5 km away from the Vertical 
Line Array (VLA) in the NW-comer of the study area were analyzed. 

Unfortunately, no sound speed data were available directly at the receiver. In order to 
define a typical sound speed profile near the VLA, several measured sound speed profiles 
closest to the receiver were averaged. The SSP used for the generation of the replica signals 
is displayed in Fig. V.l. 

Bottom depth values from the VLA towards the direction of the SUS drop, 
approximately along a 90 m isobath, were extracted from a bathymetry data base of the 
region. Bottom properties such as sound speed, sound speed gradient, density, attenuation and 
shear were not exactly known, so reasonable parameters for typical continental shelf 
sediments were used. Specifically, values chosen were compressional sound speed 
c ^ =1600m/s, no compressional sound speed gradient, a density and attenuation q = 1.5 

g /cm 3 and a^= 0.10 dB/km/Hz. Shear was neglected and the water was assumed to have a 

density of 1.0 g I cm . 

During the course of the SUS deployment, only the upper 8 hydrophones of the NW 
VLA functioned properly. The depths of these phones spanned from 30.5 m to 55 m with 
inter-element spacing of 3.5 m. The sampling frequency of the SUS data was recorded at 
3906.25 Hz. The SUS charges were 1.2 lb of TNT explosive designed to detonate at a depth 
of 18.3 m. 
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Figure V. 1 : Averaged continental shelf water sound speed profile of Mid-Atlantic Bight Field 
Study. 



As can be seen from the sound speed data displayed in Fig. V.l, there was a sound 
channel axis present in the region that spanned from about 28 m down to the bottom depth of 90 
m. The top receiver of the array was therefore in the channel while the SUS source detonation 
was above the channel. Based on the results from the synthetic data sets, this would be 
considered a very unfavourable situation, especially when much uncertainty might exist in the 
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SSP and the bottom parameters. 

Figure V.2 shows the time sequence of the recorded signal on the top hydrophone of the 
NW VLA shortly after the SUS explosion. 




Figure V.2: Real valued time-sequence of recorded SUS data 




Figure V.3: Frequency spectrum of recorded SUS data 
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The recorded signal was sampled at 3906.25 Hz over a length of 8.3886 seconds. Figure V.3 
shows the frequency spectrum of the recorded signal which illustrates that the low frequencies 
are dominant and that major contributions of the energy of the signal lie within the frequency 
band below 500 Hz. To be able to use the matching algorithms as before, the recorded signal was 
sampled and multiplied by a shifted Hanning window in order to reduce the signal to a length of 
2048 spectral points over a bandwidth of 500 Hz to 2300 Hz. 

When the search grid was drastically undersampled spatially the results were 
unsatisfactory. No localization could be made. For a spatially well sampled acoustic field and the 
shallowest receiver (30.5 m), however, the FACM algorithm made a good localization whereas 
the TACM located the source a little too shallow and a bit too far out as displayed in Figs. V.4(a) 
and (b). Attempts to improve the localization results by decreasing the minimum frequency down 
to 200 Hz rather than 500 Hz were unsuccessful. The spatial sampling of the environment proved 
to be the determining factor for a good match. Again, the FACM produced better results than the 
TACM. Attempts to localize the source with the deep receiver did not provide the desired results. 
However, for the FACM there still appears a relative peak at the correct location, and for the 
TACM the peak value appears at the correct range but on the surface instead of 18.3 m deep, as 
seen in Figs. V.4(c) and (d). The spatial sampling was the same as for the shallow receiver. 
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Figure V.4:(a) TACM, Ambiguity surface for Primer experiment, source location approximately 
at 6.48 km, 18.3 m deep, receiver at 30.5 m, spatial sampling 86X200. 
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Figure V.4 (continued) :(b) FACM. Ambiguity surface for Primer experiment, source location 
approximately at 6.48 km, 18.3 m deep, receiver at 30.5 m, spatial sampling 86X200. 
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Figure V.4 (continued): (c) FACM. Ambiguity surface for Primer experiment, source location 
approximately at 6.48 km, 18.3 m deep, receiver at 55 m, spatial sampling 86X300. 

(d) TACM. Ambiguity surface for Primer experiment, source location approximately at 6.48 km, 
18.3 m deep, receiver at 55 m, spatial sampling 86X300. 
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Finally, an attempt was made to estimate the original signal from the replica signal where 
the localization was made as a consistency check on the quality of the localization. The measured 
signal can be represented by the product A(f)g(f) where A(f) is the source spectrum and g(f) is 
the ocean transfer function, or Green's function. In the previous analysis, the replicas have been 
defined by predictions of the ocean transfer function modulated by a simple Hanning window 
over the bandwidth of interest, i.e. H(f)g’(f). In the ideal case, the predictions are exact, g'(f) = 
g(f), and a localization at the correct location would provide the true ocean transfer function 
between source and receiver. This information could then be inverted to yield a prediction of the 
source spectrum, A(f). The receiver filter characteristics were assumed to be flat (no filtering) 
over the bandwidth considered. 

In the case of the experiment considered here, the match, of course, was not perfect. 
However, the frequency spectrum of the reconstructed signal, using the method outlined above, is 
quite similar to what would be expected of a SUS source. Comparing the reconstructed spectrum 
to the spectrum of a similar 1 lb depth charge (Urick, 1975), as displayed in Figs. V.5(a) and (b), 
shows that the slopes of energy decay versus frequency match relatively well. Both drop by 
approximately 1 1 dB over one decade in bandwidth. 

The overall results of this analysis with real data are promising, especially given the huge 
uncertainty in the environment, the unfavourable relative position of source and receiver, and the 
rather unfavourable frequency spectrum of the source with the power concentrated in the lower 
frequency band. 
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Figure V.5: (a) Frequency spectrum of the reconstructed original signal, (b) Spectrum of 1 lb 
depth charge taken from Urick (Principles of Underwater Sound, 1975) 
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VI. CONCLUSIONS 



The synthetic data illustrate that the TACM is generally more useful for smaller 
bandwidths at low frequencies, which has been observed in previous work (de Kooter, 1997), 
whereas the FACM is better suited for larger bandwidths and higher frequencies. The footprint 
size for the FACM is generally larger than for the TACM, and the peak-to-sidelobe levels are 
mostly better for the FACM than for the TACM. 

Better localization results were obtained when the environmental parameters were well 
known. However, with a large enough bandwidth the FACM seems fairly robust against 
environmental mismatch. Generally an uncertainty in the SSP seems to affect the results more 
than a mismatch in the bathymetry. 

The relative position of source and receiver with respect to the SSP and bathymetry can 
affect the results significantly. Specifically, if there is a strong gradient or a sound channel, a 
rather advantageous relative positioning of source and receiver can improve the results vastly. 

Another basic observation made is that better spatial sampling of the acoustic field leads 
to improved results for both FACM and TACM. Better spatial sampling, however, increases the 
computational effort and therefore the time needed to obtain the results. An ambiguity surface for 
a spatial resolution of 3.2 m in depth and 332 m in range with an environment of size 100 m x 10 
km as used for most synthetic cases, could be created within less than 4 minutes with the 
available machine which is an operationally acceptable period of time. Certainly, with faster 
processors and more memory, this time can still be reduced. Also, the Matlab code used could at 
least partially be modified into C or C ++ which would probably further reduce the computation 
time. One could also split the environment into several spatial segments and run several routines 
at the same time, and thereby reduce the computional time needed to obtain a result and increase 
the probability of a good localization by using a better spatial sampling. A parallel run for both 
algorithms might be another way of improving results and computational time needed. 

The preliminary localization results for the Mid-Atlantic Bight Field Study are quite 
encouraging. One might next want to look at some real data with lower SNR's than for the SUS, 
which was about 8 dB in the case considered here, and evaluate the robustness of the algorithms 
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then. The results in Section IV. D show, at least for synthetic data, that even at very low SNR the 
characteristics of a signal present in noise can still be exploited for a good localization. Noise 
seems to be less critical for a good localization than the amount of environmental mismatch, one 
of the primary advantages of using autocorrelation functions. Therefore, the algorithms can be 
expected to still work well for data with a fairly low SNR. 

Another way of improving results proved to be the zero-padding of signals, especially for 
scenarios with fairly high environmental mismatch. This provides extra looks between actual 
data points and therefore provides some kind of artificially increased spatial sampling. 

Throughout this work second order moments (Chapter HI) were matched which is 
basically equivalent to the matching of instantaneous power (Hager, 1997). It was also attempted 
to match third order moments of the same quantities. The matched quantities, however, did not 
represent physically meaningful parameters and the third-order moment method was not 
successful. The same is true for the matching of just the complex pressure field data. Attempts at 
such a fully coherent matching, which represents a physically meaningful quantity, failed to 
produce good results as soon as environmental mismatch was introduced. Figure VI. 1 shows the 
poor results for Case 2, the isospeed scenario with slight uncertainty in the SSP. Another 
difficulty that occurs when using either first or third order moments is the fact that the obtained 
values for the ambiguity surface are not real but complex quantities, so that one has to compare 
either the absolute values or the phase or just real or imaginary parts for the ambiguity surface. 
Neither the matching of absolute values nor the comparison of the real parts led to successful 
results. 
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Figure VI. 1 : Case 2, Ambiguity surface for matched pressure field data, source at 3. 1 km, 
22 m deep, shallow receiver, spatial sampling 32X300. No match. 



A last point to make is that neither of the logarithmic matching methods worked well. 
Instead of enhancing later, weaker arrivals these algorithms compressed all data on the ambiguity 
surface to values between 0.9990 and 0.9999, and the footprint sizes and peak-to-sidelobe levels 
were fairly poor. As soon as environmental mismatch was introduced these methods failed 
completely. This might be due to the fact that the weaker arrivals are probably more sensitive to 
mismatches. It therefore seems that this matching method will not lead to improved results in 
possible future work considering environments similar to those considered in this work. 
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